DO ARCS REQUIRE FLAT HALO CUSPS? 
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ABSTRACT 

It was recently claimed that several galaxy clusters containing radial and tangential gravitational arcs and having a 
measured velocity-dispersion profile for the brightest cluster galaxy had to have central density profiles considerably 
flatter than those found in CDM cluster simulations. Using a simple analytic mass model, we confirm this result 
for axially symmetric mass distributions, but show that steep density profiles are well in agreement with the cluster 
requiring the flattest axially symmetric profile once even small deviations from axial symmetry are introduced. 



1. Introduction 



Are observed gravitational arcs in galaxy clusters compatible 
with the density profiles produced in CDM simulations, which 
consistently find that the dark-matter density increases towards 
| halo centres as r or s teeper? Building upon a suggestion by 
i— H iMiralda-Escudl Jl995h . lSand etalJ <2003t hereafter STSE) have 
C> recently analysed six galaxy clusters containing tangential arcs, 
CN| three of which also contain radial arcs. Apart from the lensing 
< data, the method uses constraints on the mass profile derived 
C*~) from the dynamics of the central cluster galaxies, specifically 

from its velocity dispersion profile (see also lSand et all2002l) . 
r-j . The method sets strong constraints. In a cluster showing both 
O radial and tangential arcs, the velocity dispersion measurement 
I essentially fixes the mass divided by the radius. Radial arcs con- 
O strain the slope of the projected mass profile at their location, 
<4_> and tangential arcs constrain the total mass enclosed by their ra- 
^ dial distance from the cluster centre. 

. . | Using this technique, STSE find that their sample of six clus- 
J> ters is incompatible with dark-matter density profiles propor- 
k^j tional to r~ 1 or steeper, but consistently require profiles as flat 

n nfi r -0.5 



If true, this result would be of great importance for CDM 
because such flat profiles are not found in simulations. How- 
ever, the analysis by STSE is based on the assumption of axial 
symmetry. We re-analyse their constraints here for the cluster 
Abell 383, which is the one apparently requiring the most sig- 
nificant discrepancy between simulated CDM density profiles. 
In doing so, we describe the lens model and its parameters in 
Sect. 2 and illustrate the basic reason why the method strongly 
prefers a shallow central density profile. Next, we introduce el- 
lipticity in Sect. 3 and illustrate its substantial impact. We con- 
clude in Sect. 4 that even moderate ellipticity can easily remove 
the discrepancy between the lensing observations in Abell 383 
and the typical CDM halo profiles. 

Like STSE, we use a ACDM cosmological model with matter 
density Q.q = 0.3, cosmological constant O.^ = 0.7, and Hubble 
constant Hq = 65kms _1 Mpc -1 . 

2. Lens Model 

Our lens model is adapted from STSE. It consists of a model for 
the dark-matter halo, 
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which is the profile found by Na varro et alJ i 19971) for p = 
1 . Steeper profiles with 3 — » 1 .5 are found by several othe r 
groups, among the m iMoore et all (119981) : Uin g & Sutol ^000); 
Klvp in et alJ OOOl ). The radial coordinate* = r/r s is the phys- 
ical radius divided by a scale radius r s , which STSE assume to 
be r s = 400 kpc. We adopt this value because the results are 
insensitive to it. 

In addition, the brightest cluster galaxy, assumed to be con- 
centric with the cluster, adds matter density to the c luster centre. 
Following the light profile, it is assumed to have a lJaffd (fl983) 
density profile, 

^ = W' (2) 

where now x = r/rj is the radius in units of the Jaffe radius rj. 
Fitting to the light profile of the brightest cluster galaxy, STSE 
find rj = r e /0.76, where r e is the usual effective radius of an r 1 / 4 
surface-brightness profile. For Abell 383, r e =46.75 ±2.04 kpc. 

Lensing properties are straightforwardly derived from these 
density profiles. The convergence K for the dark-matter halo is 



K h (x) =2k s jc 1_p 
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JWvithe et all 1200 II): the special ca se for the NFW profile, |3 = 
1, was derived in B arte lmannl ( 119961) . The factor K s is defined by 
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with the usual critical surface-mass density for lensing E cr . For 
the Jaffe profile, 
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CafMI!98l . with Kj = pjrjLg,. 1 . Deflection angles are derived 
from 

2 r x 

a(x) = - / yK(y)dy . (6) 
x Jo 

This needs to be computed numerically for the dark-matter pro- 
file, while 

2xacosh(x~ 1 
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for the Jaffe profile. Given the deflection-angle profile u(x), the 
radial and tangential eigenvalues are 
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Table 1 . — Fixed parameters assumed for the lensing analysis 
in this paper, taken from STSE. The two remaining parameters, 
i.e. the density profile slope p and the ratio of masses /j con- 
tributed by the dark-matter and galaxy density profiles contained 
within the Jaffe radius, are taken as free parameters. 



lens redshift 


Zd 


0.189 


source redshift 


Zs 


1.01 


effective radius 


r e 


46.75 kpc 


velocity dispersion 
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respectively (see, e.g ., ISchneider et al.l 1 19921 
Naravan & Bartelmann 1999). Radial or tangential critical 
curves are found where V = or A. t = 0. A useful relation for 
axially symmetric lens models is 

JUjc)=2k(jc)-^. (9) 
x 

In our application, source and lens redshifts together with 
the cosmological model fix the critical surface-mass density E cr . 
The Jaffe radius rj is fixed by fitting the surface-brightness pro- 
file. The mass of the brightest cluster galaxy is determined by 
the central velocity dispersion as described below. The scale ra- 
dius of the dark-matter profile is kept fixed at r s = 400 kpc. The 
only remaining parameter is the dark-matter density scale p s , 
which we express by the ratio /j between the masses contributed 
by the cluster- and galaxy density profiles within the Jaffe ra- 
dius. This ratio will effectively scale the cluster mass relative to 
the fixed galaxy mass. 

The circular velocity for an spherically symmetric density 
profile is 



For the dark-matter profile, M(r) is steeper than r, hence v rot — > 
for r — > in the dark-matter profile only. Thus, the measured cir- 
cular velocity near the centre of the brightest cluster galaxy must 
be dominated by the galaxy itself. Well within the Jaffe radius, 
the Jaffe profile is isothermal with p g °< r~ 2 . Thus, the velocity 
dispersion profile is expected to be flat close to the cluster cen- 
tre, and the velocity dispersion is approximately related to the 
circular velocity by 

Ov =V2v rot . (11) 

Through Eqs. (II Oi and dl It . the central velocity dispersion thus 
fixes the mass contained in the Jaffe profile. The assumed mass 
ratio between the masses of the dark halo and the galaxy pro- 
file, and the central slope p of the dark-matter density profile, 
thus completely determine the total density profile composed of 
a galaxy- and a dark-matter profile. 

Specifically, aiming at the cluster Abell 383 which produced 
the most significant deviation from the numerically simulated 
dark-matter profiles in the study by STSE, we adopt the fixed 
parameters listed in Tab.^ 

The convergence, deflection-angle and radial and tangential 
eigenvalue profiles k(x), a(x), |A, r (x)| and |A, t (x)| are displayed 
in Fig. [l] for p = 0.5 and mass ratio /j = 60. The figure illus- 
trates that the (projected) density profile near the cluster centre 
is dominated by the near-isothermal slope of the brightest clus- 
ter galaxy. Near the Jaffe radius, it flattens towards the central 
slope of the dark-matter density profile, and then steepens to- 
wards K °= x~ 2 as x approaches unity. Interestingly, the radial 
eigenvalue profile has two roots, indicating the presence of two 
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FIG. 1. — Radial profiles are shown for the lensing convergence K, the 
deflection angle a, and the tangential and radial eigenvalues of the lens 
mapping X t and Xr, respectively. The lens model is axially symmetric 
and composed of a central Jaffe profile and a dark matter profile with a 
central slope of — 13 = —0.5. 



radial critical curves. For Abell 383, the inner radial critical 
curve is too close to the cluster centre for having any practical 
relevance, but clusters at more favourable redshifts might show 
signatures of a double radial critical curve. This is an interest- 
ing feature of the combination of a steep, near-isothermal galaxy 
profile embedded into a relatively flat dark-matter halo. 

Having defined the axially symmetric lens model, the con- 
straints imposed by the central velocity dispersion, the radial 
and the tangential arcs are straightforwardly understood. As 
mentioned before, the central velocity dispersion is almost ex- 
clusively contributed by the mass associated with the brightest 
cluster galaxy because M/r tends to zero for r — > for the flat- 
ter dark-matter density profile. In contrast, the tangential arc is 
located at a radius which encloses a mean surface-mass density 
of unity. A first constraint thus derives from the requirement to 
have the central cluster mass dominated by the brightest cluster 
galaxy, and yet to have sufficient mass in the dark-matter halo to 
produce tangential arcs at relatively large cluster-centric radii. 
In the axially symmetric models, this is achieved by flattening 
the dark-matter density profile. 

A second constraint is imposed by the radial critical curve, 
where the derivative of the deflection angle reaches unity. For 
the Jaffe profile alone, the deflection angle is flat. It steepens 
as the total density profile becomes flatter at radii where the 
dark matter starts dominating, then flattens again as the cluster- 
centric radius approaches the scale radius of the dark-matter pro- 
file. If the dark-matter density profile is relatively flat, the in- 
crease of the deflection-angle slope occurs closer to the bright- 
est cluster galaxy than for a steeper dark-matter profile. The 
second constraint thus derives from the requirement of having a 
radial arc rather close to the brightest cluster galaxy, while the 
first constraint requires a tangential arc rather far away from the 
cluster centre. 

Figure|2]illustrates the situation. The abscissa is the mass ratio 
between dark and luminous constituents within the Jaffe radius. 
The ordinate is the cluster-centric radius in kpc. The horizontal 
bars mark the radial and tangential arc locations in Abell 383 
with their respective uncertainties. The curves in the upper and 
lower halves of the figure show the tangential and radial critical 
radii, respectively, for dark-matter profiles with four different 
central slopes, p S {0.25,0.5,0.75, 1.0}, as marked in the plot. 

A model can explain both the radial and the tangential arc if 
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FIG. 2. — Illustration of the radial and tangential critical radii for an 
axially symmetric lens model composed of a central Jaffe profile and 
a dark matter profile with four different values for the central slope, 
p = {0.25,0.5,0.75,1.0}, respectively, as indicated in the plot. The 
abscissa shows the mass ratio between the cluster component and the 
galaxy component within the Jaffe radius, the ordinate is the cluster- 
centric radius in kpc. The two horizontal stripes mark the locations of 
the radial and tangential arcs, as labelled. The curves in the upper and 
lower halves of the plot mark the tangential and radial critical radii, 
respectively. 



there is a single mass ratio for which the radial critical curve falls 
into the lower horizontal band, and the tangential critical curve 
falls into the upper horizontal band. Figure[2]confirms the result 
by STSE regarding Abell 383. Only for shallow central density 
profiles, p ~ 0.5, can the position of both the tangential and the 
radial arc be understood. For a mass ratio near 60, the lens model 
has radial and tangential critical curves in the observed ranges. 
Steeper mass profiles, e.g. the NFW model with p = 1, either 
have the tangential arc in the right range if the mass ratio is ~ 90, 
but then the radial arc is too distant from the cluster centre, or 
the radial arc location is reproduced if the mass ratio is ~ 70, 
but then the tangential arc is way too close to the cluster centre. 
The core of the problem is thus that the location of the radial 
arc requires a relatively low cluster mass, and then a steep mass 
profile forces the tangential arc too close to the cluster centre. 

3. Ellipticity 

3.1. Model and numerical results 

These conclusions are valid for axially symmetric lens models. 
As we shall show now, the situation changes considerably if de- 
viations from axial symmetry are allowed. 

As a simple model for asymmetry, we deform the lens model 
such that iso-contour lines of the lensing potential \|/ are ellipses. 
We thus introduce the radial coordinate 
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and replace \)f(x) by \\i(x). Being the gradient of \\t, the deflection 
angle a now has the components 
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where a(x) is the deflection-angle profile of the axially symmet- 
ric lens taken at x. The axially symmetric case is recovered for 
e = 0. 



Elliptical distortions of the lensing potential lead to dumbbell- 
shaped surface-mass distributions if e becomes large, e > 0.2 
day, depending on the density profile. This is certainly un- 
wanted for galaxy-sized lenses, but not necessarily for cluster 
lenses which are often highly structured. In any case, we shall 
see below that the impact of a small ellipticity e <C 1 on the lo- 
cation of the critical curves is identical for lenses with elliptical 
iso-potential curves and axially symmetric lenses embedded into 
external shear. 

Ellipticity stretches the tangential critical curves along the 
major axis of the ellipse and shrinks it along the minor axis, 
and deforms the radial critical curve in the perpendicular direc- 
tion. This implies that the cluster-centric distance of the critical 
curves now covers a range of radii. This range is surprisingly 
wide even for small ellipticities, as Fig.[3]illustrates. 
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FIG. 3. — Changes in the location of the radial and tangential critical 
curves in response to a moderate elliptical distortion of the lensing po- 
tential. The figure is arranged to resemble Fig. [3] but showing results 
for p = 0.5 and (3 = 1 .0 only for clarity. Three curves are shown for 
each profile slope. The central curve shows the critical radius for axial 
symmetry, the other two curves indicate the range of radii of the critical 
curves for an ellipticity parameter of e = 0. 1. Even moderate ellipticity 
widens the radial ranges such that the observed radial and tangential 
critical radii can quite easily be reached even for steep profiles. 

The figure is arranged in the same way as Fig. |2] and also 
specialised for the cluster Abell 383. The radius from the clus- 
ter centre is plotted against the mass ratio between cluster and 
brightest cluster galaxy for different central slopes p of the clus- 
ter density profile. The observed locations of the tangential and 
radial arcs are marked as horizontal bars. For clarity, we now 
show results for two values of p only, p 6 {0.5, 1}, but illus- 
trate with three curves for each P the radial range covered by 
the tangential and radial critical curves for a small ellipticity pa- 
rameter e = 0. 1 . For a fixed mass ratio, the radial range is given 
by the vertical distance between the top and bottom curves of 
the same type. Conversely, for a fixed radius, the horizontal dis- 
tance between the left and right curves of the same type show 
the range of mass ratios for which arcs at that radius can be pro- 
duced somewhere along the respective critical curves. The cen- 
tral curves reproduce the critical radii for the axially symmetric 
case. 

Evidently, the impact of the small ellipticity e on the location 
of the critical curves is quite large. For the shallow profile, P = 
0.5, and a mass ratio of 60, the tangential critical radius ranges 
from 27 to 70 kpc, while the axially symmetric result is 45. For 
the same mass ratio and density-profile slope of p = 0.5, the 
radial critical radius ranges between 3 and 19 kpc around the 



axially-symmetric value of 7 kpc. Thus, even the low ellipticity 
of e = 0. 1 makes the location of the tangential and radial critical 
curves vary by about a factor of 1 .5 and more than a factor of 
two, respectively, around the cluster. Interestingly, the radial 
ranges for tangential and radial critical curves now overlap even 
for the steep profile with p = 1 . For a mass ratio of ~ 80, the 
observed positions of both the radial and the tangential arcs fall 
within the ranges allowed by the model. 

3.2. Analytic description 

These results can be reproduced analytically. For coordinate 
axes aligned with the elliptical iso-potential contours, the mini- 
mum and maximum values of the critical radii occur on the axes. 
The intersection of the tangential critical curve with the xi-axis 
satisfies 

a(x) 

1- - v J — = 0. (14) 
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Expanding around the solution for the critical radius in the axi- 
ally symmetric case, which satisfies 
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FIG. 4. — Relative shift per unit ellipticity, Sac/ (ex), of the radial and 
tangential critical curves caused by external shear or internal ellipticity. 
Results for two values of the central density-profile slope (3 are shown. 
For example, for a mass ratio of 60, (3 = 0.5 and £ = 0.1, the radial 
critical curve shifts by ~ 100%, the tangential curve by ~ 40%. 



and assuming e <C 1, we find the amount 8x by which the tangen- 
tial critical curve is shifted on the xj-axis relative to the axially 
symmetric critical radius. Repeating the calculation for the jc 2 
axis, we obtain 

ex 

where the + and — signs apply to the x\ and xi axes, respec- 
tively, and Kt is the convergence at the tangential critical radius 
of the axially symmetric model. Repeating this analysis for the 
radial critical curve yields 



8x r = ± 



ex 



2(K r -l+*£c) ' 
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where K r and kJ. are the convergence and its radial derivative at 
the radial critical radius of the unperturbed lens model. Figure^ 
shows the relative shifts for unit ellipticity, i.e. 8x t r /(ex), for the 
two central slopes p s {0.5, 1}. 

The figure shows that, for a mass ratio of 60 and a central 
density slope of p = 0.5, the relative shifts for e = 0.1 are of 
order 100% for the radial critical curve, and 40% for the tangen- 
tial critical curve, confirming the numerical results illustrated in 
Fig-El For the steeper profile with p = 1, the respective relative 
changes are of order 40% and 20% for the same mass ratio and 
ellipticity. 

It is straightforward to show that the first-order results d!6t 
and dl7> remain valid if the lens model is not itself deformed, 
but embedded into an external shear y. In that case, y simply 
replaces e in these equations which are otherwise unchanged. 

Equations dl6> and illi show that the amount by which el- 
lipticity or shear shift the critical curves depends highly sensi- 
tively on the slope of the convergence profile k(x). For a sin- 
gular isothermal sphere, for instance, K = 0.5 at the location of 
the tangential critical curve, which is at x = 1 . Thus, 8 t = ±e 
in this case. Flatter profiles, however, have K closer to unity in 
Eq. dl6l because the tangential critical radius encloses a mean 
convergence of unity. The relative shift of the tangential critical 
curve is thus amplified for flatter density profiles, as illustrated 
in Fig.0]for the two choices of p. As the figure also shows, the 
situation is more extreme for the radial critical curve. Flatter 
profiles thus much more sensitive to external shear or internal 
ellipticity than steeper profiles. 



3.3. Numerical Example 

We illustrate the analytic results obtained above using a nu- 
merically simulated galaxy cluster as a lens. It is located 
at redshift z = 0.3 and has a mass of M = 5 x IOU/j^Mr. 
The cluster was taken from a large-scale numerical simulation 
of the ACDM model with parameters Q.q = 0.3, £l\ = 0.7, 
Ho = 70kms~ 1 Mpc -1 and normalisation Og = 0.9. It is one 
of the clusters produced by the GIF collaboration tKald eich. 
1999 1) used in earlier related stu dies (iBartelmann et al.ll l998: 
Meneghe tti et all l2000. 2003b a). The particle mass in the simu- 
lation is 1 .4 x 10 10 Ii~ 1 Mq and the gravitational softening length 
was set to 30 hr kpc. 

We first use the cluster as it is, i.e. with the asymmetry 
and substructure produced by the simulation, and then progres- 
sively smooth and circularise it by computing its azimuthally- 
averaged density profile, subtracting it from the cluster, smooth- 
ing the residual density distribution by a varying amount, and 
finally adding the axially symmetric density profile back on the 
smoothed residuals. Using a normalised smoothing kernel, this 
procedure preserves the total mass and the mean density profile 
of the cluster. 

Figure[5]illustrates the effect of the progressive smoothing on 
the critical curves. The ragged lines show the critical curve of 
the original clusters, while the approximately elliptical, smooth 
curves are the critical curves of the cluster after smoothing. Evi- 
dently, smoothing makes the critical curves shrink considerably, 
and, more importantly for our discussion, the radial range cov- 
ered by the critical curves narrows substantially while the mean 
density profile remains entirely unchanged by construction. Ar- 
guments based on the cluster-centric distance of radial and tan- 
gential arcs thus need to take the detailed cluster structure into 
account. 



4. Summary 

Galaxy clusters containing radial and tangential gravitational 
arcs and a brightest cluster galaxy with a measured velocity dis- 
persion apparently require significantly flatter density profiles 
than obtained in numerical simulations of CDM cosmologies. 
As shown here, this is essentially caused by the large observed 
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FIG. 5. — Critical curves of a simulated cluster in various stages of 
smoothing. As described in the text, the smoothing procedure con- 
serves the total cluster mass and its density profile by construction. The 
ragged line is the original critical curve. Even moderate smoothing 
makes the critical curves shrink and considerably reduces the radial 
range where tangential and radial arcs can be found. 



cluster-centric distances of tangential arcs, which require fairly 
flat density profiles given the central constraints of radial arcs 
and the velocity-dispersion measurement. Using a simple ana- 
lytic mass model, we can confirm the results by STSE, provided 
the tensing mass distribution is axially symmetric. 

Allowing deviations from axial symmetry, the results radi- 
cally change. We have chosen to introduce asymmetry by ellip- 
tically distorting the lensing potential, but showed that embed- 
ding the axially symmetric lens into external shear has identical 
consequences in the limit of small ellipticity or shear. Based on 
these results, we have shown that the particular cluster which 
most significantly required a flat density profile in the analysis 
by STSE, Abell 383, is well compatible with an NFW profile 
((3 = 1) even for the small ellipticity of e = 0. 1 . 

Critical curves caused by flat density profiles are e xtremely 
sensitive to distortions, as demonstrated in Sect. l3.2.1 and illus- 
trated in Fig. |4] Shifting tangential critical curves by 20% to 
40%, and radial critical curves by 50% to 100% even with an 
ellipticity or shear of only 0.1 is possible in particular for the 
profiles as flat as advertised by STSE. This is also the reason 
why the analysis of cluster ellipticity carried out by STSE them- 
selves concluded that ellipticity had a negligible effect on their 
results: their lens model used components with isothermal den- 
sity profiles which are much less sensitive to external shear or 
distortions, as illustrated by Eqs. ( II 61 and ( II 71 . 

We conclude that radial and tangential arcs in clusters do not 
rule out central density profiles as steep as found in CDM sim- 
ulations once effects of asymmetry and shear are taken into ac- 
count. We will extend our analysis towards numerically simu- 
lated clusters in a forthcoming paper. 
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